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THEORETICAL CALCULATION OF THE EFFECT OF THE FUSELAGE 
ON THE SPANWISE LIFT DISTRIBUTION ON A WING 
By Martin Zlotnick and Franklin W. DIederich 

SUMMARY 


A method is presented for calculating the effect, of the fuselage 
on the spanwlse lift distribution on a wing by an application of a 
conformal -mapping procedure to the simplified lifting-surface theory. 

This method is applicable to any symmetrical wing-fuselage configuration. 

At low angles of attack the method is shown to approach the theoretically 
correct limit for wing-fuselage configurations having a wing of vanishingly 
small aspect ratio and for those having a wing of infinite aspect ratio. 
Some calculations made by this method are compared with experimental 
data. A computing scheme is given in which the calculation for the infi- 
nitely long fuselage and the correction for the finite fuselage length 
are made separately. A numerical example is given in an appendix to 
illustrate the entire procedure. 


INTRODUCTION 


A knowledge of the aerodynamic loading on the wing of an airplane 
is required for the stress analysis of the wing and for the prediction 
of the flying qualities of the airplane. Experimental determination of 
this loading is tedious and is resorted to in very few cases. Satisfac- 
tory theoretical methods for the calculation of the spanwise lift distri- 
bution on the wing alone have been developed, but the problem of predicting 
the effect of the fuselage on the spanwise lift distribution on the wing 
has not yet been satisfactorily treated. Experience has indicated that 
the effect of the .fuselage tends to_be smallest low angle dof “attack and, 
for”"the lack of ahv Tbetter informati on u sual-lT y^neglecte d. However, 
refinements in the methods of stress analysis place a premium on the 
accurate prediction of the aerodynamic loading applied to the wing and, 
hence, on the effect of the fuselage on that loading. Also, for swept 
wings, indications are that the effects of the fuselage on the spanwise 
lift distribution of the wing are not as small as for unswept wings. 
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Consequently, development of a theoretical method for predicting these 
effects is very desirable at this time. 

In est ima t ing the effect of the fuselage on the spanvise lift distri 
bution on a wing, it is necessary in general to consider the following: 

(1) The fuselage causes a change In the longitudinal velocity in the 
vicinity of the wing. This effect Is hereinafter referred to as the 
"inflow effect." 

(2) If. the fuselage is at an angle of attack relative to the free 
stream, it changes the flow about the wing in planes normal to the free 
stream. This effect is hereinafter referred to as the "cross flow effect 

(3) The fuselage has a blocking effect^ since there can be no flow 
normal to Its surface. 

The inflow effect, which is a result of the finite length of the 
fuselage, changes the longitudinal velocity in the vicinity of the wing 
and, hence, the lift on the wing. This change Is not large for a slender 
fuselage but may be important when the fuselage Is relatively thick, so 
that a substantial change Is caused in the local longitudi n a l velocity. 
For the same reason this effect may be important for a wing -nacelle 
combination. 

The cross-flow caused by the fuselage changes the component of free- 
stream velocity normal to the fuselage axis and also affects the dowawash 
flow produced by the wing. The blocking effect of the fuselage Is always 
present even If the fuselage Is an Infinite cylinder alined with the free 
stream, in which case the other two effects are not-^present . 

•» 

In the following analysis the blocking and cross-flow effects will 
be treated by an application of a conformal -mapping procedure to the sim- 
plified lifting-surface theory, and the inflow effect will be treated as 
a separate correction. The conformal -mapping procedure for calculating 
the blocking and cross-flow effects is derived in a somewhat intuitive 
fashion from Multhopp f s method (reference 1) and is not the oretic ally 
Justified because it applies a conformal transformation to a flow which 
Is not strictly two-dimensional. However, the method presented herein 
(referred to as the new method) is shown to be in good agreement with 
experimental results and to give the same result as more exact theory for 
two limiting cases - namely, for wing -fuselage- configurations having a 
wing of infinite aspect- ratio and for those having a wing of vanishingly 
small aspect ratio. 



NA.CA RM L51J19 


3 


SYMBOLS 
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chord 

mean chord (s/b) 

section lift coefficient 

section lift -curve slope (per radian) 

geometric angle of attack of wing 
geometric angle of attack of fuselage 

section lift; also length of body 
mass density of air 
free -stream, velocity 

longitudinal component of free-stream velocity (V U) 
local longitudinal velocity 

local, velocity increment - ij 

circulation 
lateral coordinate 
longitudinal coordinate 
vertical coordinate 
wing area 
wing span 

aspect ratio (b 2 /s) 
radius of body 
diameter of body 
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A 

X 

v 3cA 

w n 

Wn f 


7 

7 


c* 


angle of sweep 
taper ratio 

downwash velocity at three-quarter-chord line 

component of velocity normal to plane of wing 

component of velocity normal to plane of wing due to 
fuselage 

dimensionless circulation (4-r/hV) 
dimensionle ss circulation ( 4r/Ev ) 

chord made dimensionless with respect - to transformed semi- 

— (v;) 


complex- coordinate in real plane (z +-iy) 
complex coordinate in transformed plane (z += iy) 

real part of conformal-transformation factor 

Superscript: 

* indicates coordinate made dimensionless with respect to 

semi span 



Subscripts: 

v refers to wing 

wf refers to wing of the wing-fuselage combination 

A bar indicates transformed coordinate^— or quantity made dimension- 
less with respect to transformed coordinate.- .. 
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DESCR IPTION OF THE METHOD 
Analysis for the Infinitely Long Fuselage 


In order to calculate the fuselage cross-flow and blocking effects 
on the spanwise lift distribution on a wing, a conformal -mapping pro- 
cedure is applied herein to the Weissinger L-method (reference 2) in the 
same manner that Multhopp applies a conformal -mapping procedure to the 
lifting -line method. 

Inasmuch as the concept of conformal mapping is based on two- 
dimensional considerations and inasmuch as the simplified lifting-surface 
method cannot, in general, be reduced to a two-dimensional problem as can 
lifting-line theory, the method outlined herein is not rigorously justi- 
fiable. However, for wings cf very low aspect ratio, this method is 
found to be justifiable because in those cases the simplified lifting- 
surface method can be reduced to a two-dimensional problem. The results 
obtainable by the method under consideration will also be shown to con- 
stitute a very close approximation to the theoretically correct results 
in the case of wings of infini te aspect ratio. Therefore, there is reason 
to believe that this method will furnish useful answers for wings of inter- 
mediate aspect ratios. 

The integrodifferential equation of the simplified lifting-surface 
method is given in reference 2. In the presence of the fuselage the 
normal -flow component on the wing w n is 


w n = Va + w^ 


Cl) 


where 


<x Is the local angle of attack on the wing and 


v. 




Is the Incre- 


mental normal flow on the wing caused by the fuselage. The lift distri- 
bution is now calculated for a wing-fuselage configuration obtained from 
the given one by a conformal transformation (in planes normal to the free- 
stream velocity) in which the cross section of the fuselage is a vertical 
slitj so that there is no fuselage Interference. The velocity components 
normal to the free -stream velocity and parallel to the plane of symmetry 
are then multiplied by the conformal -transformation factors and the sim- 
plified lifting-surface equation is written for the wing in physical space^ 
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but in terms of the transformed coordinates. Thus, in physical space, 
the normal -flow component-caused by the fuselage is 


= Va f ■ Va f (y) 


( 2 ) 


where (y) is the angle of attack of the fuselage} also, the downwash 
velocity at the three -quarter-chord line becomes 


w 3cA = 



( 3 ) 


where 



in Multhopp'B 
the preceding 


lift in terms 


is the real part' of the conformal-transformation factor as 

treatment of the wing -fuselage combination. As a result of 
relations the integrodifferential equation expressing the 
of the wing and fuselage geometry becomes 


a(y*) 


+ a^y* ) r(§|) - 1 


p 1 drj* 

W ^ J.i y* - n* 

— "- 1 " i 


+ 
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»(?*) - <*(?») . 1 r 1 aL. as + 

_/d£\ ^ ^ J-i d Ti* y* - n* 

K «y 


1. 


8ir 


2 



(5) 


Comparison of equation (5) with equation (3*0 of reference 2 shows 
that the procedure for calculating the lift distribution on the wing- 
fuselage combination is only a slightly modified form of that required 
for the wing alone. 


Computing Procedure for a Midwing Configuration 

In order to illustrate the method, the case of a cylindrical fuselage 
with a midwing is outlined in this section and a numerical example for a 
configuration with a "shoulder high" wing is given in the appendix. Other 
fuselage cross sections and wing locations may be treated in the same 
manner provided the correct transformation factor is used. (See 
reference 1 . ) 

The lateral and vertical coordinates in physical space, y and z, 
respectively, are combined into one complex coordinate 


t = z + iy 


( 6 ) 


Similarly, these coordinates in the space obtained by transforming each 
plane in physical space normal to the free-stream velocity are combined 
into a complex coordinate 


t = 2 + iy 




(7) 


For a cylindrical fuselage, the relation between the two complex coordinates 
is v.. 

, 0 - 

1 A 1 *'* 

\ ~ 






SA- 


w~'> 




■■U 


% *' 


v.\ ^ r 


£ = 5 + y 


6 

*<» 


I 

it 


A ^ 




( 8 ) 




* 



8 


NACA RM L51J19 


where a is the radius of the fuselage. The calculation for the span- 
wise lift distribution is then made as follows: 

(l) In the transformed plane, the coordinate y for a midwing con- 
figuration is 


y = : 

so that when y = a, y = 0. When 

1=1 

b 

Also, since the longitudinal coordinate- is not - changed by the transformation, 

x = x (9b) 



(2) Since 


& = 1 _ S? 

£2 


( 10 ) 


the real part of the conformal -transformation. factor is 



(n) 


for z = 0, the- -vertical, position of the wing. 

(3) If the wing is unswept;, equation (5) is solved in the usual manner 
for the transformed wing, except that- a - is divided by 


The- 
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difference a - ctf is 0 when the wing is not set at an angle of inci- 
dence with the fuselage axis . 

(4) If the wing is swept, another method must he used for treating 
the transformed wing, since the transformed swept wing has a curved 
quarter-chord line (see fig. 1 ), a case for which Weissinger's method 
cannot he applied. 

The simplest procedure in this case is prohahly the use of horseshoe 
vortices (about 10 centered on the quarter-chord line on each semispan 
appear to he a satisfactory number) . By calculating the downwash induced 
hy each vortex at 10 points on the three -quarter-chord line of the semi- 
span and equating the sum to the slope of the wing at the three -quarter- 
chord line, a set of 10 simultaneous equations Is obtained for the 10 
unkn own strengths of the vortices. (See the appendix.) In calculating 
the downwash due to the vortices, the tables in reference 3 are very 
helpful . 

(5) The relation between the dimensionless circulation on the wings 
in physical and transformed space is 


7 



( 12 ) 


( 6 ) The dimensionless circulation 7 is plotted at the lateral 
coordinate in the physical plane which, from step ( 2 ), is 


y* - 



(13) 


where 


y* 


= -X- 
b /2 


If , for instance, the lift distribution is known on a rectangular 
wing alone, the lift distribution for the combination of this wing and 
a fuselage can be obtained quite readily. The results of the calculation 
for the wing alone can be assumed to be the results for the transformed 
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wing except for a small-correction which can he made by applying lifting- 
line theory. This procedure may be used if the ratio of the span to the 
fuselage diameter is large-so that - b/b (equation ( 9 a)) is nearly unity. 
The transformation is then made back into the physical plane so that the 
lift is plotted as indicated in step ( 6 ). If 7 ^ is the dimensionless 

circulation on the wing-fuselage combination and 7 ^ Is the dimension- 
less circulation on the wing alone, then 


7vf(y*) - 


— llWi 

7 w (y*)l — — 
l w y a + 


A + 2 


( 1 *) 


The correction applied in equation (l4) is based on the lifting -line 
solution for the elliptical wing 


c l 


2 it 


A 

A + 2 


a 


(15) 


and the identity which holds in the case of a rectangular wing 


7 



(16) 


Correction for the Finite Length of the Fuselage 

In representing the fuselage by an Infinite cylinder the inflow 
effect Is neglected. The lateral distribution of the Induced velocity 
is shown in figure 2 for two fairly typical bodies} the body with the 
slenderness ratio of 3 approximates a nacelle and the body with the slender- 
ness ratio of 7 approximates a slender fuselage. An increase of the effec- 
tive free-stream dynamic pressure is the Imme diate result of the induced 
longitudinal velocity. If the increase in the velocity Is small and does 
not have a large lateral gradient (as, for example, in the case of the 
slender body In fig. 2) and the wing has a high aspect ratio, the effect 
on the spanwise lift distribution may be calculated as follows: 

If the local velocity at a given spanwise station is U(l + 6 ) where 
8 « 1, the lift in the region of increased velocity Is 

I = CjC^pU 2 ^ +-8) 2 (17) 
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or 


Z ^cgpl^fl + 2B) 


( 18 ) 


but the lift in the uniform flow l Q is 


Zq — CjC^p 


-U 2 


( 19 ) 


therefore. 


Z * (1 + 25) Z Q 


( 20 ) 


Conse quen,tly, the lift calculated by the method outlined in the preceding ' 
section can be corrected approximately for the inflow effect by multiplying 
it by the factor 1 + 28 appropriate to the given fuselage and any span- 
wise station. 

It should be noted that the value of 5 on the surface of a given 
fuselage varies -with the longitudinal position of the wing and also with 
the vertical location of the wing when the fuselageis at an angle of 
attack. Prom figure 10 of reference 4 , more lift is seen to be expected 
on a high-wing than on a low-wing configuration if the wing is located 
ahead of the position of the fuselage maximum diameter. It should also 
be noted that, since the wing is of finite thickness, the longitudinal 
velocity induced by the fuselage may be different on the wing upper sur- 
face from that on the lower surface} this difference tends to increase or 
decrease the circulation on the wing and, hence, to change its lift 
accordingly. The inflow correction given before may therefore be applied 
only when the wing section is thin relative to the fuselage diameter. 
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DISCUSSION OF METHOD AND CALCULATED RESULTS 
Comparison of the Results of the New Method for the Limiting Case of 
Vanishingly Small Aspect Ratio with the Results of 
Linearized Low-Aspect -Ratio Theory 


In order to obtain the spanwise lift distribution on the wing- 
fuselage combination for the limiting case of vanishingly small aspect 
ratio. It Is only necessary to multiply the value of 7 (y*), the loading 
coefficient for the wing alone, by the factor b/b and to plot the 
result 7 ^ at - the points on the wing-fuselage in the physical plane 

which correspond to the points for the wing alone in the transformed 
plane. (See steps (5) and ( 6 ) of the section entitled "Computing Pro- 
cedure for a Midwing Configuration.") The relation between the lateral 
coordinates in the physical plane and the lateral coordinates in the 
transformed plane is given in equation ( 9 ). The solution for the wing 
alone of vanishingly small aspect ratio may be obtained from reference 5 
as 

7 W (! f*) = 4 1/1 - y * 2 (21) 


which can be used for the wing-fuselage combination by setting 


?vf (y*) = | ?w(y*) 


( 22 ) 


Figure 3 shows a comparison of the results obtained by using the 
new method for the limiting case of vanishingly small aspect ratio with 
the results of Spreiter's method (reference 6 ) which Is theoretically 
correct for the limiting case . As would be expected, no discrepancy Is 
apparent because, in this limiting case, the problem of calculating the 
lift distribution is essentially a two-dimensional one (in planes per- 
pendicular to the stream) j thus, the procedure of transforming the physical 
space plane, for plane into a distorted space Is rigorously Justifiable. • 

It should be noted that - the results shown In figure 3 apply only to 
configuration witn wings of very low aspect ratio (A < 1) where the flow 
is essenti all y two-dimensional in planes normal to the free stream and, 
therefore, the trend of decreased lift- on the wing in the presence of 
the fuselage would not necessarily apply to configurations with wings of 
moderate or high aspect-ratio. 
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Comparison of the Results of the New Method at the Limiting Case of 
Infinite Aspect Ratio with the Results Obtained from 
Considerations of an Idealized Case 

The calculation made by Lennertz (reference 7} to determine the span- 
wise lift distribution on a combination of a constant -strength vortex 
with a cylinder alined with the free stream can be modified slightly to 
give the lift distribution when the cylinder axis is placed at an angle 
of attack with the free stream. In this case a vertical velocity com- 
ponent w n normal to the axis of the cylinder must be considered in 

addition to the longitudinal velocity U considered by Lennertz so that 
the resultant stream velocity V is 



+ 


U 2 


(23) 


and the angle of attack of the cylinder is 


, _i w n 
a, , = tan — 

cyl *“* u 


(210 


In a potential flow this additional upward velocity cannot change the 
lift on the cylinder; however, this velocity can give rise to an additional 
lift on the part of the wing (represented by the vortex) adjacent to the 
fuselage . 

The magnitude of this increase in the lift will be calculated by 
treating the components of the lift Zy and l v , due, respectively, to 

the components of the free-stream velocity U and w n , as follows: 

The component of the lift Zy on the wing which acts normal to the 
cylinder axis is 

Zu = pUT (25) 

The component of the lift Z w which acts parallel to the cylinder axis 
and opposite in sense to the velocity U Is 



kv / 1 U- 


V /°or 
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■where w n ^l + is the velocity of the vertical flow at the vortex 

line in the presence of the cylinder. But, as a result of equation (24), 


Z v = pur tan 


4 • 0 


The component of 2rr normal to the free -stream velocity is 


‘•U = i U 


and the component of l v normal to the free-stream velocity is 


l v 1 = ,l v sin a 
w n n 


The resultant lift l normal to the free-stream velocity is therefore 


l = Ijj cos a + 2 W ^ sin a 


= pur cos a + pur — — (1 + 

K COS a\ 


J*Ul ♦ sin a 
cos a\ y2 t 


l a 2 

= pVr 1 + — sin^a 


( 30 ) 
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since U = V cos a. Now, if the angle of attack a is small, the term 

a 2 0 

— sin a is negligible j therefore 

y 2 


l 5 ; PVT (31) 


which is the result for the wing alone. Consequently, as a first-order 
approximation, the lift distribution on the unswept wing of infinite 
aspect ratio is unaffected by the presence of the fuselage. The same 
result is given by equation (l4) for the limiting case of the wing- 
fuselage combination having a wing of infinite aspect ratio. 


Comparison of the Calculations with Experiment 

The results of calculations made by the method presented in this 
paper are compared in figures 4 and 5 with experimental results obtained 
from references 8 , 9, 10, and from unpublished data. In figure 6 the 
results of a calculation by Multhopp*s method are shown and compared with 
experimental data from references 8 and 9- For the unswept wing (figs. 4 

and 6 ), the fuselage was assumed to be a Rankine solid ^ = 7.0 (fig. 2) 

and for the swept wing the fuselage was assumed to be an ellipsoid of 

revolution — = 10 . 0 . 
h 


Experimental and theoretical results shown in figures 4 and 5 indi- 
cate that the swept wing is affected more by the fuselage attached to it 
than is the unswept wing, although the fuselage radius of the swept wing 

b b 

configuration is only 0 .l|: as compared with a fuselage radius of 0 . 19 — 

on the unswept wing configuration. The theoretically calculated increase 
in lift is attributed largely to the decreased sweep at the root of the 
transformed wing (fig. 1 ). 

It should be noted that, in calculating the inflow effect, the wing 
thickness was assumed to be negligible in each case, although the wing 
thickness was 10 percent of the fuselage diameter and 15 percent of the 
fuselage diameter for the unswept and swept configurations, respectively. 
This assumption leads to a value for the lift which is slightly less than ■ 
the correct value. 

Figures 4 and 6 show comparisons of experimental results of refer- 
ences 8 and 9 with results of calculations made by the method of this 
paper and the method of Multhopp. A large discrepancy exists at the 
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root of the wing between ■ Multhopp * s results and the results of experiment, 
whereas the agreement of experimental results with the results calculated 
by the method of this paper is fairly good along the entire span. The 
large lift calculated by Multhopp at~the root of the wing is attributed 
by him to his neglect of the effect of- the fuselage on the dovnwash of 
the bound vortex, since he calculates the downwash infinitely far behind 
the wing rather than at the three -quarter-chord line. 


CONCLUDING REMARKS 


A method has been presented for calculating the lift distributions 
in incompressible flow over swept— or unswept wings with fuselages. This 
method is more generally applicable than methods presented heretofore 
and the calculations made by this method are in fair agreement with 
experimental results. Although further comparison with experimental lift 
distributions is desirable before the- method is applied generally, one 
conclusion may be drawn from the available results: The presence of a 

slender fuselage does not have an important effect on the lift distri- 
bution on an unswept wing of moderate aspect ratio, but a larger change 
in the lift distribution on a wing in the presence of a fuselage may be 
anticipated if the wing is swept - ; 


Langley Aeronautical Laboratory 

National Advisory Committee for* Aeronautics 
Langley Field, Va. 
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APPENDIX 

ILLUSTRATIVE EXAMPLE OF THE CALCULATION OF THE SPANWISE LIFT 
DISTRIBUTION ON THE WING IN THE PRESENCE OF THE FUSELAGE 
Geometry of the Wing-Fuselage Combination 


A plan view and elevation of the wing-fuselage combination con- 
sidered in the illustrative example are shown in figure 7- Some of 
the pertinent geometrical parameters are: 


A = 8.02 a = 0.1C& 
, 2 

\ = 0.45 k = 10 

b - 

A = 450 z F = 0.05§ 


Calculation of Coordinates of Transformed Quarter-Chord Line 

The coordinates of the wing are transformed according to the 
formulas (see equation (8) of text and reference l) 


y = y (1 - 


y 2 


V 


fi- 

b/2 b 


where y = y = and £ = 1 - 


*2 


„2 „2 

y + z F * 

_*2 


1 + 


.2- 


At a given value of 
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The coordinates of the quarter-chord lines of the transformed and 
un transformed wing are as follows: 


y* 

y 

h /2 

** 

II 

<'|xi 

ro 1 

O.O 8667 (wing root") 

0 

0 

.1 

.020 

.0133 

.2 

.154 

.1133 

.3 

.270 

.2133 

.4 

• 379 

.3133 

• 5 

.484 

.4133 

.6 

.588 

.5133 

• 7 

.692 

.6133 

.8 

.795 

.7133 

.9 

.898 

.8133 

1.0 

1.000 

.9133 


Calculation of Chords on the Transformed Wing 

The formula for the chord length at a given spanwise station may Toe 
written 


c 


* 


4 

A(l + X) 


[l - (1 - X)y3 


therefore 


=£- 4 c* =4r [l 

h/2 h h A(l + X) L 



For the plan form "being treated 


c = 1.01 X 4 / 

h/2 8.02 X 1.45 \ 


o.55y*) 
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The chords may he tabulated as follows: 


y 

b/2 

y* 

c 

b/2 

0 

0.0866 

0.331 

.05 

.120 

.325 

.15 

.198 

.310 

*25 

.283 

*293 

.35 

•374 

.276 

.45 

.468 

.258 

.55 

.561 

.240 

.65 

.660 

.221 

• 75 

.758 

.203 

.85 

.854 

.184 

.95 

.951 

.166 


Calculation of Lift Distribution on the Transformed Wing 

The lift on the wing is calculated by replacing the wing by a set 
of 20 lifting horseshoe vortices located on the quarter-chord line at 
the spanwise stations y* =±0.05, ±0.15, ±0.25, ±0-35, ±0.45, ±0.55, 
±0.65, ±0.75, ±0.85, and ±0.95* The downwash induced by these vortices 
on the three-quarter-chord line of the wing at these same spanwise 
stations is equated to the angle of attack of the wing at the three- 
quarter-chord line. This procedure yields a set of 10 simultaneous equa- 
tions of the form ( reference 3 ) : 


03c/lj.( x *>y*) 




The numbers are arranged so that 



(1 < n < 10) 


and 


yn* = -y*n- 


n-10 


( 11 ^ n ^ 20) 
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The longitudinal location of 1 the quarter-chord line of- the wing is 

is the 


represented by x n * which corresponds to y n * and 
loading coefficient of station x n *,yn*- 


c Jla) 

J >/2 In 


Since tlie load is symmetrical, 


CM -(£*) 

\b/2 Jn \b/2/ n -10 


(11 % n £ 20) 


and the equations may"' therefore he written 


“30/^** , y *) = 45 - 6o ^ (l^) m Gm t X * " ( y * “ ym *)] (A1) 


where 




(1 ^ n ^ 10) 


and 


G m = F n + F n+ i 0 (l - n = 10) 


The values of F may he most conveniently obtained from the tables in 
reference 3 . In the notation of reference 3* 


Ax v =20 (x* - x m *) 
Ayv = 20 (y* - y m *) 


Equation (Al) may be written for the 10 different values of y*. 
Since the values of a and F are known in each of the 10 equations, 

cc l 

the equations may be solved simultaneously for the values of — — 

b/2 


CC Z 

The values of — — 
b/2 

and are multiplied by T>/ b 


so obtained are given in the following table 
CC Z 

to give — : 

b/2 
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_z_ 

b/2 

y* 

cc 7 
L _o, 

b/2 

cc 7 

V2 

0.05 

0.120 

1.085 

1.074 

.15 

.198 

1.087 

1.077 

.25 

.283 

1.078 

1.067 

-35 

•374 

1.046 

1.036 

.^5 

.468 

1.000 

.99 0 

• 55 

.561 

.947 

.938 

.65 

.660 

• 895 

.886 

• 75 

.758 

.832 

.824 

.85 

.854 

.752 

.744 

• 95 

.951 

.622 

.616 


Calculation of the Inflow Effect 

The value of u/U on the surface of a body of revolution is 
obtained very simply from equation (22) of reference 4. For an ellips- 
oid of revolution, this equation reduces to 


(!f ♦ v] 

where 

v - Vx : (i ) 2 

and 

<fcCv> - I ^Se^-rr) - 1 

When Z/h is 10 and x Is 0, then v = 1.0Q5 and 

Ql = 0*5025 log e §^§§| - 1 = 2.01 

and 

u = 1^05 (0 * 01 x 2,01 + 1 - 005} = 1 - 02 
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But, since ij = 1 + 8, then 8 must- "be 0.02 on the surface of the body. 

Since 8 is small, the lateral variation of 8 can be taken to be 
the same as the lateral variation of 8 in figure 2 without any sig- 
nificant loss in accuracy in the final result. 

As shown in the text, the lift at a section is multiplied by the 
value of 1 + 28 which corresponds to the spanwlse location of the 
section to correct for the inflow effect. The values of 1+28 and 
of the corrected lift distribution in the case- of this example are as 
follows: 


* 

y 

1 + 28 

cc ? 

L a 

b/2 

( corrected for 
inflow effect) 

0.120 

1.04 

1.117 

.198 

1 . 04 — 

1.120 

.283 

1.04 

1.110 

.374 

I.03 

1.067 

.46 8 

1.03 

1.020 

.561 

1.02 

.956 

.660 

1.02 

.904 

.758 

1.01 . 

.832 

.854 

1.01 

• 752 

• 951 

1.00 

.6 16 


This lift distribution is plotted in. figure 5 
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Dimensionless lateral coordinate, y/a 


Figure 2,- Lateral velocity distribution for two Rankine solids 
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Wing alone 


f Wing with fuselage x. 
(Theory of reference 6) ^ 


O Theory of this paper 


Dimensionless lateral coordinate, y " 

Figure 3. - Theoretical lift distribution on a combination of a fuselage 
with a wing of low aspect ratio, a* = 0.3. 


Wing with fuselage- 


Wing alone 


□ Experiment, wing with fuselage (reference 9) 
O Experiment, wing alone (reference 8) 


Dimensionless lateral coordinate, y 
Figure It.- Comparison of theoretical and experimental lift distributions 


on unswept wing with and. without fuselage . 
a* = 0 . 19 ; a ■» 1*°. 


h.$3 X = 1# A 
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Figure Comparison of theoretical and experimental lift distribution 
on a swept wing with and without fuselage . . A = 83 X = A = h$°; 

a* = 0 . 1 ; a ■ U. 7 °. 



Fig-are 6.- Calculation for wing with and without fuselage. by Multhopp's 
method. A = U.5>> X = lj A = 0°; a* = 0.2; a' = I 4 0 . 





Figure 7 


- Plan view and elevation of wing-fuselage configuration 
considered in illustrative example. 
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